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ORIGINAL ARTICLE 



Differentially expressed genes and signalling pathways 
are involved in mouse osteoblast- 1 ike MC3T3-E1 cells 
exposed to 17-p estradiol 

Zhen-Zhen Shang\ Xin Li^, Hui-Qiang Sun^'^, Guo-Ning Xiao\ Cun-Wei Wang^ and Qi Gong^ 

Oestrogen is essential for maintaining bone mass, and it has been demonstrated to induce osteoblast proliferation and bone formation. 
In this study, complementary DNA (cDNA) microarrays were used to identify and study the expression of novel genes that may be 
involved in MC3T3-E1 cells' response to 17-p estradiol. MC3T3-E1 cells were inoculated in minimum essential media alpha (a-MEM) 
cell culture supplemented with 17-p estradiol at different concentrations and for different time periods. MC3T3-E1 cells treated with 
10~^ mol L~^ 17-p estradiol for 5 days exhibited the highest proliferation and alkaline phosphatase (ALP) activity; thus, this group was 
chosen for microarray analysis. The harvested RNA was used for microarray hybridisation and subsequent real-time reverse 
transcription polymerase chain reaction (RT-PCR) to validate the expression levels for selected genes. The microarray results were 
analysed using both functional and pathway analysis. In this study, microarray analysis detected 5 403 differentially expressed genes, 
of which 1 996 genes were upregulated and 3 407 genes were downregulated, 1 553 different functional classifications were identified 
by gene ontology (GO) analysis and 53 different pathways were involved based on pathway analysis. Among the differentially expressed 
genes, a portion not previously reported to be associated with the osteoblast response to oestrogen was identified. These findings clearly 
demonstrate that the expression of genes related to osteoblast proliferation, cell differentiation, collagens and transforming growth 
factor beta (TGF-p)-related cytokines increases, while the expression of genes related to apoptosis and osteoclast differentiation 
decreases, following the exposure of MC3T3-E1 cells to a-MEM supplemented with 17-p estradiol. Microarray analysis with functional 
gene classification is critical for a complete understanding of complementary intracellular processes. This microarray analysis provides 
large-scale gene expression data that require further confirmatory studies. 
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INTRODUCTION 

Oestrogen, a major sex steroid hormone, has garnered considerable 
attention because of its importance in bone mass maintenance and the 
efficacy of hormone therapy in combating osteoporosis. 'Bone mass is 
maintained by balancing osteoblast-mediated bone formation and 
osteoclast-mediated bone resorption. Oestrogen deficiency leads to 
an increased rate of bone remodelling and tilts the balance between 
bone resorption and formation in favour of the former, resulting in 
bone loss and osteoporosis.^"^ Oestrogen deficiency is not only the 
dominant pathogenic risk factor for postmenopausal osteoporosis, 
but also the key factor that contributes to continuous bone loss in 
aging men.'' 

It has been well established by both in vivo and in vitro studies that 
oestrogen inhibits osteoclastic bone resorption.^ This phenomenon 
that has been variously attributed to the suppression of cytokine 
production in the bone microenvironment, which includes interleu- 
kin-1, interleukin-6, macrophage colony-stimulating factor, receptor 
activator of nuclear factor kappa-B ligand (RANKL) and tumour 



necrosis factor-ot,^"^ and to the induction of apoptosis in osteoclasts.* 
Bone remodelling is a complex ceU-ceU interaction; osteoblasts and 
osteoclasts are regulated by one another. The receptor activator of the 
nuclear factor kappa-B (RANK)-RANKL-osteoprotegrin (OPG) 
pathway was an important discovery that explained the mechanism 
by which the osteoblast lineage influences the activity and differenti- 
ation of osteoclasts.' RANK is expressed on the surface of osteoclast 
precursors, RANKL is found on the surface of osteoblasts, and OPG, a 
'decoy receptor' molecule, is released by osteoblasts. The activation 
and differentiation of osteoclast precursors into mature (active) osteo- 
clasts requires the binding of RANKL to RANK. The RANK-RANKL 
interaction is inhibited by OPG, which competitively binds to 
RANKL.' Previous studies have revealed that 17- (3 estradiol increases 
osteoprotegrin mRNA and protein levels in human osteoblastic cells 
in a dose- and time-dependent manner.'" Alternatively, a paracrine 
mechanism was found in which oestrogen affects osteoclast survival 
via upregulation of the Fas ligand in osteoblasts, leading to the apop- 
tosis of osteoclast precursors." 



^Department of Prosthodontics, College of Stomatology, Shandong University, Jinan, China; ^Department of Stomatology, The Affiliated Hospital of Medical College Qingdao 
University, Qingdao, China and ^Shandong Provincial Key Laboratory of Oral Biomedicine, Shandong University, Jinan, China 

Correspondence: Dr HQ Sun, Department of Prosthodontics, College of Stomatology, Shandong University, Wenhuaxi Road 44-l#, Jinan 250012, China 
E-mail: whitedove69@163.com 
Accepted 18 December 2013 



Responsive genes and pathways to 17-|t estradiol in MC3T3-E1 cells 

ZZ Shang etal 

143 



Oestrogen's ability to prevent bone loss is also attributed to its ability 
to stimulate bone formation.^''^ Oestrogen receptors (ERs) on osteo- 
blasts have been identified as potential target cells for oestrogen 
replacement therapy/^''' Oestrogen helps preserve bone mass by 
increasing osteoblast proliferation and function. Postmenopausal 
women exposed to relatively high doses of oestrogen demonstrate 
sustained stimulation of osteoblast function.'^ Furthermore, oestrogen 
prevents apoptosis of osteoblasts via the repression of apoptotic gene 
expression, thereby extending the life span of osteoblasts/^ Oestrogen 
deficiency also downregulates the transcription of the gene for insulin- 
like grovrth factor 1, reducing bone formation as a result.'^ Various 
intermediary pathways have been found that participate in the oestro- 
gen-mediated biological responses of osteoblasts, such as transforming 
growth factor (TGF)-P, the Wnt signalling pathways and apoptosis. 

Although oestrogen's various mechanisms of action on bone cells 
have been widely investigated, the studies are primarily focused on one 
or a small number of key factors; studies providing a complete picture 
from high-throughput data and signalling pathways via gene chips are 
rather limited. In this study, an in vitro experiment combined with a 
gene microarray technique was performed to identify novel genes that 
may be involved in MC3T3-E1 cells' response to oestrogen. This large 
dataset can be examined for bioinformatics approaches that allow the 
examination of biological systems rather than alterations in individual 
genes, which wiU facilitate further understanding of the molecular 
mechanisms of oestrogen that affect osteoblasts. 

MATERIALS AND METHODS 

Cell culture 

MC3T3-E1 osteoblast-hke cells (ATCC, Manassas, VA, USA) were 
cultured in minimum essential media alpha (a-MEM) supplemented 
with 10% foetal calf serum (Invitrogen, New York, NY, USA) and 1% 
penicillin/streptomycin (Invitrogen, New York, NY, USA) at 37 °C in 
humidified air with 5% CO2. The cell culture medium was changed 
every 2-3 days. 

Cell viability 

Cell viability was determined using the 3-(4,5-dimethylthiazol-2-yl)- 
2,5-diphenyltetrazolium bromide (MTT) assay kit (Sigma-Aldrich, St 
Louis, MO, USA) according to the manufacturer's protocol. The cells 
were seeded in each well of a 96-well plate with 200 |iL of culture 
media. The cells were incubated in the presence of 17- (5 estradiol at 
different concentrations (10"'°, 10""*, 10"** and 10"^ mol-L"'), and 
each concentration of 17-P estradiol was maintained for 1, 3, 5 and 7 
days. A group that did not receive oestrogen treatment served as the 
control. At the end of the incubation period, the cells were incubated 
for 4 h with the MTT reagent at a final concentration of 0.5 g-L~' 
followed by washing with phosphate-buffered saline (pH 7.4) and the 
addition of dimethyl sulfoxide. Complete dissolution was achieved 
after gentle shaking. The absorbances of the resulting solutions were 
recorded at 490 nm using a microplate spectrophotometer. Each con- 
dition was performed in quintuplicate. 

Alkaline phosphatase assay 

Alkaline phosphatase (ALP) activity was determined via enzymatic 
assay. Experiments were conducted 24 h after cells were seeded onto 
a six-well chamber. MC3T3-E1 osteoblast-like cells were treated with 
17-P estradiol at different concentrations (10"'", 10~', 10~* and 
10~^ mol-L~'), and each concentration of 17-P estradiol was main- 
tained for 3, 5 and 7 days. A group that did not receive oestrogen 
treatment served as the control. After treatment, the cells were rinsed 



with phosphate-buffered saline and then lysed in buffer containing 
10 mmol-L"' Tris-HCl, pH 7.4, 0.2% Triton X-100 (Beijing Native 
Biological Technology, Beijing, China). The cell lysates were centri- 
fuged, and the supernatants were used for the assays according to the 
manufacturer's protocols. The optical density of p-nitrophenol at 
520 nm was determined spectrophotometrically, and the results were 
normalized to the protein content of the sample, which was quantified 
using the bicinchoninic acid (BCA) protein assay kit (Keygenbio, 
Nanjing, China) in accordance with the manufacturer's protocol. 

Extraction of total RNA 

The most appropriate concentration and incubation time of 17-P 
estradiol based on the results of the ALP assay and the MTT assay were 
applied to another group of MC3T3-E1 cells, which were cultured and 
processed as previously described. A Trizol Reagent kit (Invitrogen, 
New York, NY, USA) was used to extract total cellular RNA. The 
extracted RNA was stored at -70 C The purity of the isolated RNAs 
was determined by measuring the optical density (OD) value (A260/ 
A280) using a spectrophotometer (Agilent, Shanghai, China). 

Genome-wide oligonucleotide microarray analysis 

Microarray analysis was performed using a Mus musculus 12X135k 
v2 array, which consists of 135 000 amino acid-modified 60-mer 
oligonucleotide probes representing 44 170 genes (Kangchen Bio 
Shanghai, China). After preparing 1 |ig of DNase-treated total RNA, 
fluorescent dye (Cy3-dCTP)-labelled cDNA produced via RNA amp- 
lification and subsequent enzymatic reaction were then hybridized 
into an array using the NimbleGen Hybridisation System 4 (Roche 
NimbleGen, Madison, WI, USA). Finally, the arrays were scanned 
using a GenePix 4000B single-channel scanner (Molecular Devices, 
Sunnyvale, CA, USA), and the data were extracted from the obtained 
images using GenePix Pro v6.0 (Molecular Devices, Sunnyvale, CA, 
USA). GeneSpring vll.O software (Agilent Technologies, Santa Clara, 
CA, USA) was used to analyse the genes that were differentially 
expressed between the experimental and control groups. The genes 
that were consistently altered in both arrays with differences in mean 
expression ratios that were greater than twofold on average were 
selected as differentially expressed genes. The microarray results were 
analysed via both gene ontology (GO) and pathway analysis using 
NimbleScan v2.5 Software (Roche NimbleGen, Madison, WI, USA). 
Each condition was performed in triplicate. 

Quantitative real-time reverse transcription polymerase chain 
reaction 

To confirm the microarray results, three representative genes were 
chosen at random for analysis with quantitative real-time reverse tran- 
scription chain reaction (RT-PCR). cDNA was prepared from 2 |ig 
DNase-treated total RNA from either the test or the control sample 
using the First Strand Superscript II Kit (Invitrogen, New York, NY, 
USA). Quantitative RT-PCRs were performed using the DNA Master 
SYBR Green I Kit and a LightCycler (Roche Diagnostics, Mannheim, 
Germany) according to the manufacturer's protocols, and the results 
were analysed using LightCyler software version 3.5 (Roche Diagnostics, 
Mannheim, Germany). Glyceraldehyde-3-phosphate dehydrogenase 
gene {Gapdh) was used as a housekeeping gene. AH data were reported 
as values normalized to Gapdh. A mathematical model reported by 
Pfaffl was employed to analyse the relative expression ratios of these 
genes.'^ The relative expression ratio was determined using the formula 
^cpi-cp2)^^^^^^cp3-cp4)^^^^^^ in which E represents the quantitative 

RT-PCR efficiency and CP is the crossing point. In the same samples. 
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Table 1 Sequences of primers used in quantitative real-time RT-PCR 


Gene symbol 


Encoded protein 


Forward primer (5'-3') 


Reverse primer {5'-3') 


Gapdh 
Lefl 
Vangl2 
Rasgrpl 


Glyceraldetiyde-3-ptiosptiate dehydrogenase 
Lymphoid enhancer binding factor 1 
van Gogh-like 2 

Ras guanine-releasing protein 1 


GTTGTCTCCTGCGACTTCA 
CGGCGTTGGACAGATCAC 
CGGGCTATTCCTACAAGTCG 
AACCACTTACCTGAAGCCCAC 


GCCCCTCCTGTTATTATGG 
GGGTAGGGTTGCCTGAATC 
CCGTGACCACCGTTGTTG 
TCTTGCACTCGAACACCACC 


RT-PCR, reverse transcription polymerase chain reaction. 



Gapdh was continuously expressed across the different physiologic time 
points. Primer sequences for each gene target are shown in Table 1. 
After reverse transcription, relative quantification of the selected genes 
was performed using a PGR Thermal Cycler (TaKaRa Biotechnology, 
Dalian, China). The reaction mixtures were incubated in the Applied 
Biosystems Rotor-Gene 3000 Real-time PGR System (Gorbett Research, 
Sydney, Australia). The Peltier thermal cycler experimental protocol 
was performed as follows: denaturing program (95 d for 5 min), amp- 
lification and quantification program with 40 cycles (95 tZ, 10 s; 60 "C, 
15 s; 72 "C, 20 s; 80 "C, 15 s with a single fluorescence detection), melting 
curve program (72 0-99 "G, with a heating rate of 0.2 tZ-s"' and 
continuous fluorescence detection). Three parallel reactions were per- 
formed for each gene. 

Additional microarray information 

This microarray study followed the minimum criteria provided by the 
microarray experiment guidelines.^" Detailed protocols for RNA 
isolation, amplification, labelling and hybridisation can be provided 
by the authors upon request. 

Statistical analysis 

All experiments were repeated at least three times using different 
specimens, and the representative data of ALP was presented as the 
means ± standard deviation. Significant differences among the MTT 
and ALP results were analysed using various analyses. StatisticaUy 
significant differences in gene expression between each pair of experi- 
ments and the corresponding control samples were analysed using 
Lowess Normalisation and a paired f-test. Greater than 2.0-fold higher 
and less than 0.5-fold lower alterations in expression with a P<0.05 
were considered significant. 

RESULTS 

Cell viability based on the MTT assay 

There were no significant differences in cell viabflity between the 
experimental groups and the control group on the first day 
(P>0.05). After 3 days, the cell viability of the groups treated with 
lO"** and 10"** mol-L"' 17-P estradiol was higher than that of the 
other two groups (10"^mol-L"^ and 10"'" mol-L"') (P<0.05); there 
was no significant difference between the 10~^ mol-L~ ' group and the 
10"'" mol-L~' group (P>0.05); After 5 days, the group treated with 
10~^ mol-L~' 17-P estradiol exhibited significantly higher cell viabi- 
lity than the other three experimental groups, and all four experi- 
mental groups performed significantly better than the control group. 
After 7 days, aU groups exhibited trends that were identical to those 
that were detected after three days. When we compared the different 
time points of treatment with 10"'* mol-L~' 17-P estradiol, we found 
that the group that received 5 days of 17-P estradiol treatment exhibi- 
ted the most significant difference in cell viabflity compared with the 
control group, which was significantly different compared with the 
other time points and concentrations (Figure 1). 



ALP activity assay 

As shown in Table 2, the group treated with 10~^ mol-L~' 17-P estra- 
diol for 5 days displayed the highest level of ALP activity. There were 
no significant differences between the experimental groups and the 
control group on the third day (P>0.05). After 5 days of treatment, the 
ALP activity of the 10~^ mol-L~' 17-P estradiol group exhibited sig- 
nificantly higher ALP activity compared with the other three experi- 
mental groups. After 7 days of treatment, there were significant 
differences between the experimental groups and the control group 
(P<0.05), but there were no significant differences detected among 
each of the experimental groups (P>0.05). 

RNA quality assurance via spectrophotometry 

The OD A260M280 ratio was approximately 2.0, which demonstrated 
that the RNA extracted from the ceUs was not degraded. 

Gene microarray analysis 

In this study, microarray analysis revealed 5 403 differentiaUy 
expressed genes in the oestrogen-treated group compared with the 
control group. Of the differentially expressed genes, 1 996 were upre- 
gulated and 3 407 genes were downregulated. The differentiaUy 
expressed genes related to bone metabolism in the oestrogen-treated 
group compared with those of the control group are presented in 
Table 3. Apoptosis and osteoclast differentiation-related genes were 
downregulated. Genes responsible for osteoblast proliferation and 
differentiation were upregulated. The genes that encode coUagens 



■1d iSd 5d iZd 




Control 10-10 io-9 10-8 iq-^ 



Figure 1 Graph of MTT results for each group at different oestrogen concen- 
trations and treatment periods (Mean±standard deviation, n=20). The MTT 

levels of the control group and the experimental groups at different oestrogen 
concentrations (10" 10"^, 10"**, 10"' mol-L"' 17-p estradiol) and different 
durations. Error bars denote standard deviation. *P<0.05, compared with the 
control group; ^P<0.05 compared with the other experimental groups of the 
same duration, but different oestrogen concentrations. MTT, 3-(4,5- 
dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide. 
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Table 2 ALP results for each group at different oestrogen concentrations and treatment periods 



/(U- 



Time 


xm-^° mol-L"i 


XIO"^ mol-L"' 


XlO"**mol-L"' 


XlO"' mol-r' 


Control 


Day 3 
Day 5 
Day 7 


0.013 1 ±0.00 6 
0.104 5±0.01 5 
0.030 7±0.00 8* 


0.014 8±0.00 6 
0.162 8±0.01 1* 
0.033 1 ±0.00 8* 


0.016 6±0.00 4 
0.276 5±0.01 1*" 
0.033 7±0.00 6* 


0.015 6±0.00 3 
0.227 5±0.01 1* 
0.022 1±0.00 4* 


0.012 9±0.00 4 
0.086 2±0.00 6 
0.026 7±0.00 1 



ALP, alkaline phosphatase; s.d., standard deviation. 
Mean±standard deviation, n=15. 
*P<0.05 compared with the control group. 

" P<0.05 compared with the experimental groups of the same duration, but different oestrogen concentrations. 

Table 3 Changes in gene expression in the oestrogen-treated group compared with the control group 



GenBank accession ID 



Gene product 



Gene symbol 



FC absolute 



Pvalue 



145 



Upregulated 

Cell proliferation 



14184 


Fibroblast growth factor receptor 3 


FgfrS 


16.611 1 


0.000 022 67 


Cell differentiation 










16779 


Lamlnin beta 2 


Lamb2 


3.945 2 


0.019 6 


ECM-receptor interaction 










12814 


Collagen type XI, alphal 


Colllal 


2.002 4 


0.006 8 


1280 


Collagen type II, alphal 


Col2al 


3.985 9 


0.024 1 


12826 


Collagen type IV, alphal 


Col4al 


2.755 4 


0.003 5 


12832 


Collagen type V, alpha2 


Col5a2 


6.640 1 


0.012 2 


Gap junction 










432530 


Adenylate cyclase 1 


Adcyl 


2.200 4 


0.012 9 


Wnt signalling pathway 










93840 


van Gogh-like 2 


Vangl2 


15.860 0 


0.001 8 


TGF-p signalling pathway 










12162 


Bone morphogenetic protein7 


Bmp7 


7.647 3 


0.036 4 


Downregulated 










Osteoclast differentiation 










20846 


Signal transducer and activator of transcription 1 


Statl 


2.432 0 


0.000 4 


Apoptosis 










26420 


Mitogen-activated protein kinase 9 


Mapk9 


8.415 0 


0.000 1 


MARK signalling pathway 










19419 


RAS guanyl releasing protein 1 


Rasgrpl 


30.970 0 


0.000 9 


Wnt signalling pathway 










56637 


Glycogen synthase kinase 3 beta 


GskSb 


11.830 0 


0.000 8 


Matrix metallopeptidases 










17390 


Matrix metallopeptidase-2 


Mmp2 


2.455 8 


0.031 9 


17387 


Matrix metallopeptidase-14 


MrnpM 


3.449 5 


0.009 8 


ECM, extracellular cell matrix; MARK, mitogen-activated protein kinase; TGF-p, transforming growth factor beta. 






Table 4 Principal biological processes based on GO analysis 








GO term 


Gene symbols 


Fold enrichment 


Pvalue 


Regulation 


Cell differentiation 


Eng, Lamb2, Tiel, Rps6ka2 


2.624 9 


7.963 2X10"=^' 


Up 


Cell proliferation 


Fgfr3, Calca. Cavl, Cdhl3 


2.506 2 


3.671 3X10"^^ 


Up 


Intracellular signal transduction 


Dok4. Cav4, Mapk8ip2, Ndstl 


2.478 8 


3.202 6X10"^'^ 


Up 


Cell communication 


Mapllc3a, Fgfrl, Spryl, Sphkl 


2.441 8 


1.833 2X10"^^ 


Up 


Regulation of gene expression 


Bcl6b, Cdknlc, Hoxa2, FoxpS 


1.658 1 


3.964 3X10"" 


Up 


Response to hormone stimulus 


Thra, Bcarl, Pik3rl, FoxoB 


3.226 8 


2.465 3X10"<^ 


Up 


Programmed cell death 


Bcl2lll, Atp7a, Jun, Mapk9, Ticaml 


2.024 6 


2.485 6X10"^^ 


Down 


Apoptosis 


BCI2I11, Atp7a, Jun, SehnbS, Mapk9 


2.029 6 


3.020 9X10"^^ 


Down 


kp kinase/NF-Kp cascade 


Iraki, Tlr4, Maltl, Tmed4 


2.001 2 


1.238 4X10"'^ 


Down 



GO, gene ontology. 



(types II, IV, V, XI) and matrix metallopeptidases (MMPs) 14 and 
MMP2 were also upregulated. 

Among the 1 553 different functional classifications related to GO 
analysis, a typical fiinctional group is presented in Table 4, including 
GO terms, gene members, fold enrichment and the P value. Upre- 
gulated GO categories included cell differentiation, cell proliferation, 



intracellular signal transduction, cell communication and response to 
hormone stimulus. The downregulated categories included programmed 
cell death, apoptosis, and the IkP kinase/NF-KP cascade. 

Pathway analysis indicated that 53 pathways were upregulated and 
67 pathways were downregulated. Upregulated pathways included the 
Wnt signalling pathway, the extracellular cell matrix (ECM) receptor 
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Table 5 Pathway analysis of differentially expressed genes 


Pathways 


Genes 


Fisher Pvalue 


Regulation 


Cell adhesion molecules 


Cd2, CdhS, Esam 


5.341X10"' 


Up 


Wnt signalling pathway 


Vangl2, Axin2, Btrc, Camk2b, Camk2d 


1.567X10"' 


Up 


Focal adhesion 


Bcarl, Cavl, Colllal, Collla2, Col2al 


4.703X10"= 


Up 


ECM-receptor interaction 


Colllal, Col2a2, Col4al, Col5a2, Col6al, Hspg2, Sdc3, Agrn 


6.074X10"' 


Up 


Calcium signalling pathway 


Adcyl, Atp2a3, Cacnalc, Calm3 


0.001 1 


Up 


Gap junction 


Adcyl, AdcyS, Gjal, Lparl, Pdgfa 


0.004 2 


Up 


MARK signalling pathway 


Cacnalc, Cacna2dl, Cacnbl, CacnbS 


0.007 3 


Up 


Osteoclast differentiation 


Ifnar2, Chuk, Csflr, Cyba, Cyld, Fos 


2.021X10"= 


Down 


Wnt signalling pathway 


Gsk3b, Camk2a, Ccndl, Ctnnbl, Ctnnbipl, Fbxwll 


0.018 1 


Down 


MARK signalling pathway 


Rasgrpl, Cacnala, Chuk, Duspl, DusplG, Dusp4 


0.038 5 


Down 


Adherens junction 


Ctnnbl, Iqgapl, Lefl, Nik, Rac2, Smad2 


0.047 1 


Down 


Cell cycle 


Anapcl, Ccndl, Ccne2, CdcMa 


0.011 1 


Down 


ECM, extracellular cell matrix; MARK, mitogen-activated protein kinase. 






Table 6 Results of RT-PCR analysis of selected differentially expressed genes 








Fold chanf 




Pvalue 


GenBank accession ID 


Gene name Gene symbol Gene chip 


RT-RCR Gene chip 


RT-RCR 



16842 Lymphoid enhancer factorl Lefl -11.77* -2.510 0.002 6 0.001 6 

93840 van Gogh-like 2 Vangl2 15.86* 10.50 0.000 6 0.008 9 

19419 Ras guanyl releasing protein 1 Rasgrpl -30.97* -35.85 0.000 9 0.000 2 

RT-RCR, reverse transcription polymerase chain reaction. 
* P<0.05. 

'-' stands for downregulated (the oestrogen-treated group relative to the control group). 



interaction, the gap junction and the mitogen-activated protein kinase 
(MAPK) signalling pathways. There was some overlap between the 
downregulated and upregulated pathways, including the Wnt signalling 
pathway and the MAPK signaUing pathway. Additionally, osteoclast 
differentiation and adherens junction pathways were downregulated. 
The corresponding pathways are presented in Table 5. 

Confirmation of differential expression of selected genes via real- 
time RT-PCR analysis 

The microarray results were validated using real-time PGR assays of 
selected genes. The findings confirmed the results obtained from the 
microarray analysis, despite slight disparity. Collectively, the results of 
the quantitative PCR demonstrated the reliability of the microarray 
analysis (Table 6). 

DISCUSSION 

After menopause, reductions in the level of circulating oestrogen result 
in bone loss. Previous studies have confirmed that oral maxillofacial 
bone loss is closely associated with systemic osteoporosis. Oestrogen 
deficiency increases variability in the tissue mineral density of the 
alveolar bone surrounding the teeth during early remodelling " and 
decreases the osseointegration index in the implant/bone interface.^^ 
In addition, 17-P estradiol has been shown to increase the osseointe- 
gration index at the implant/bone interface in osteoporotic rats.^^ 

17-P estradiol is a type of artificially synthesized exogenous oestro- 
gen with an activity and function similar to that of endogenous oes- 
trogen, and it is often used to study the function of oestrogen. The 
MC3T3-E1 cell line, established from newborn mouse calvaria, is a 
useful model system for studying the mechanism of osteogenesis, as it 
possesses ERs (ERa and EPP) and has been reported to retain the 
capacity to differentiate into osteoblasts.^'' In this study, we applied 
different concentrations and durations of 17-P estradiol to MC3T3-E1 



cells. Viability and proliferation were evaluated using the MTT assay, 
which may be a more sensitive method for assessing osteoblast pro- 
liferation because it measures cell viability via the determination of 
mitochondrial dehydrogenase activity.^^ The proliferation of MC3T3- 
El cells reached its peak when subjected to 10~^ mol-L~' 17-P estra- 
diol for 5 days. As the treatment period was extended further, the 
proliferation activity of the MC3T3-E1 cells tended to decline. 
These results are consistent with the experimental evidence reported 
previously. ^''"^^ The differentiation of osteoblasts was assessed via ALP 
activity, an early marker of osteoblast differentiation.^^ After 3 days of 
treatment, the ALP activity of MC3T3-E1 cells was low because the 
cells were still in the proliferation stage. The ALP activity of MC3T3- 
El cells reached its peak when the cells were subjected to 10~^mol-L~' 
17-P estradiol for 5 days; as the treatment period was extended, the 
ALP activity of MC3T3-E1 cells tended to dechne. The MTT and ALP 
results indicate that the response of the MC3T3-E1 cells to 17-P estra- 
diol varies according to the concentration and duration of oestrogen 
and the cell stage. Thus, the 10~* mol-L~' 17-P estradiol group was 
chosen for the microarray analysis. 

Genome-wide analysis via microarrays is a useful and complemen- 
tary approach to genetic studies for the identification of genes and 
corresponding pathways that contribute to a given phenotype. An ana- 
lysis of the genes vnth differential expression in MC3T3-E1 cells in 
response to 17-P estradiol could yield information about potential 
signal transduction pathways and candidate genes. Biological process 
analyses enable the further understanding of the differentially expressed 
genes. In this study, a total of 40 genes were involved in hormonal 
stimulus responses, and 174 genes were involved in intracellular signal 
transduction. Additionally, a number of routine biological processes 
exhibited elevated expression, including those responsible for cellular 
proliferation {Fgfr3, Calca, Apoe, Cavl, Eng, Vashl and CdhlS) and cell 
differentiation {Lamb2, M112, Rps6ka2, Cavl and Heyl). Fibroblast 
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growth factor receptor 3 (FGFR3) is an important regulator of bone 
development and remodeUing. Studies have indicated that FGFR3~'~ 
mutations in mice result in decreased cortical bone thickness and osteo- 
malacia.^' Therefore, these processes may imply an underlying molecu- 
lar mechanism by which 17-P estradiol promotes the prohferation and 
differentiation of MC3T3-E1 cells. 

Two hundred and six genes involved in apoptosis {Bcl2lll, Atp7a, 
Jun, Mapk9, Serpinb9, Tepl, etc.) were downregulated based on the 
biological processes analysis. Bradford et al?" induced apoptosis of 
G-292 human osteoblasts and this activity was inhibited by pre-treat- 
ment with estradiol. They concluded that oestrogen combats osteo- 
porosis by preventing osteoblast apoptosis and that a key event in this 
process is the repression of apoptotic gene expression. The downre- 
gulation of apoptosis-related genes may provide another functional 
effect of 17-P estradiol on osteoblasts. 

It was reported that the beneficial effects of oestrogen on osteoge- 
nesis are caused in part by oestrogen's ability to suppress osteoclasto- 
genic cytokine production by osteoblasts.^^ Thirty genes {Stall, Ifnar, 
Ifngr, Chuk, Csflr, Mapk9, Ncfl) related to osteoclast differentiation 
were also found to be downregulated in this experiment. Signal trans- 
ducer and activator of transcription 1 (STATl) is a critical regulator of 
both osteoclastogenesis and osteoblast differentiation in skeletal frac- 
ture healing. Studies reveal that STATl negatively regulates osteoblast 
differentiation by inhibiting the nuclear translocation of RUNX2 and 
by suppressing osterrx transcript levels and promoter activity.^'"^^ 
These results indicate that STATl interferes with osteoblast-mediated 
bone formation. Interferons (IFNs) are important factors in the 
immune system. Studies of bone destruction associated with rheu- 
matoid arthritis have highlighted the importance of the interaction 
between the immune and skeletal systems. Previous studies have indi- 
cated that IFN-P modifies human osteoblast function by inhibiting 
ECM synthesis, eventually resulting in delayed bone formation and 
mineralisation in the early phase of osteoblast differentiation.^"* 
Oestrogen represses osteoclast activity via the downregulation of genes 
related to osteoclast differentiation, which maybe another mechanism 
by which oestrogen promotes osteoblast proliferation. 

The RANK-RANKL-OPG system is an essential signalling pathway 
involved in bone ceU-ceU communication, and it plays important 
roles in bone remodeUing. OPG, a key modulator of the function of 
the RANK ligand, plays a significant role in the physiological regu- 
lation of and pathological alterations in bone metabolism.^'' It has 
been demonstrated that OPG increases bone mass when the recom- 
binant protein is administered to normal rats,^ "' and the deletion of the 
mouse Tnfrsfllb results in severe osteoporosis resulting from increased 
osteoclast activity.'^ In this study, the expression of Tnfrsfllb was not 
significantly altered. Previous studies have indicated that the protein 
and mRNA expression of osteoprotegrin were significantly increased 
after 24 h of estradiol treatment, while smaller increases were observed 
with estradiol compared with untreated primary human osteoblast 
cells at extended time points.^'' The lack of a significant difference in 
the expression level of Tnfrsfl 1 b may be attributed to the cell system, or 
it may be caused by the extended duration of oestrogen treatment in 
our study. 

The Wnt signalling pathway is one of the principal pathways par- 
ticipating in osteoblast differentiation and bone formation. Studies 
have revealed that Wnt pathway activation enhances osteoblast and 
osteocyte survival in vifro,^^"'"' and the Wnt pathway was active in 
bone regeneration sites.'" Oestrogen has also been implicated in the 
upregulation of Wnt protein and Frizzled receptor expression,''^ and it 
induces the differentiation of uncommitted osteoblast progenitors by 



stimulating both Wnt and bone morphogenetic protein (BMP) signal- 
ling in a kinase- dependent manner in vitro and in vivo^^ In this study, 
microarray analysis identified the differentially expressed genes involved 
in Wnt; the expression of synthase kinase 3 beta (Gsksp) reduced by 10 
times, and the expression of van Gogh-like 2 [Vangll) increased more 
than 10-fold. These results were confirmed via RT-PCR. GSK3P is a 
negative regulator of canonical Wnt signalling that phosphorylates the 
serine and threonine residues of P-catenin, leading to the degradation of 
P-catenin via a proteasome-dependent mechanism.'''' Wnt signalling 
disrupts the GSK3p-Axin-APC-P-catenin complex via poorly under- 
stood mechanisms that may involve the recruitment of GBP/frequenfly 
rearranged in advanced T-cell lymphoma (Prat). GBP/Frat displaces 
Axin from GSK3P, leading to P-catenin release and stabilisation. 
Gaur et al.'^^ found that the canonical Wnt pathway promotes bone 
formation via the P-catenin/T-ceU factor (TCF)-mediated activation of 
the master osteogenic transcription factor Runx2, which drives 
mesenchymal cells to an osteogenic hneage. In addition to the canonical 
Wnt signalling pathway, Vangll, a non-canonical Wnt signalling gene, 
was also affected by the oestrogen treatment of MC3T3-E1 cells. Alysia 
found that a Vangll mutation in the loop-tail mouse uterine epithelium 
displayed altered cell polarity"*^ concomitant with changes in cytoskele- 
tal actin and scribble locahsation. 

In addition to the Wnt signalling pathway, the MARK pathway plays a 
key role in extracellular signal transduction to initiate cellular responses. 
MARK pathways can relay, amplify and integrate signals from a diverse 
range of stimuU and elicit appropriate physiological responses, including 
cellular proliferation, differentiation, development, inflammatory res- 
ponses and apoptosis. Studies have demonstrated that the MARK sig- 
nalling pathways activate and phosphorylate the osteoblast-specific 
transcription factor Cbfal in MC3T3-E1 cells and play an important 
role in the regulation of osteoblast-specific gene expression."*^ Other 
studies suggest an important role for the extracellular signal-regulated 
kinase (ERK)-MAPK pathway in the regulation of osteoblast differenti- 
ation and foetal bone development.''^ Alternatively, 17-P estradiol was 
found to induce osteoblast proliferation via the G protein-coupled 
receptor GPR-30, which is partly mediated by MARK activation. 

Connexin 43 (Cx43), which is expressed in virtually all types of bone 
cells, is the major component of gap junctions and plays important roles 
in signal transduction. Previous studies have revealed that the dynamic 
balance of bone mass was broken by reduced Cx43 expression rates in 
the osteoblasts and increased Cx43 expression rates in the osteoclasts of 
ovariectomized rats.''' Reportedly, oestrogen upregulates Cx43 express- 
ion in the human myometrium.''" Moreover, Cx43 expression was 
shown to be downregulated by oestrogen deficiency in myocardial cells 
in ovariectomized rats.^' Recent studies indicate that oestrogen not only 
increases Cx43 expression and the function of gap junction intercellular 
communications, but also enhances the mechanosensitivity of MLO-Y4 
cells to mechanical loads via ERs.''^ In this study, the expression level of 
Gjal increased, which is consistent with previous studies. These findings 
all point to a possible Unk between oestrogen, Cx43 and Cx43-based gap 
junction intercellular communication. An elucidation of the interaction 
between oestrogen and Cx43 will deepen our understanding of the 
effects of oestrogen on bone metabolism. 

BMPs are the most potent regulators of osteoblast differentiation 
among the local factors. Osteogenic protein 1, also referred to as BMP- 
7, belongs to the TGF-P superfamily. BMP-7-deficient mice exhibit 
Polydactyly and occasional abnormalities of the ribs.^"* Studies indi- 
cate that BMP-7 induces the expression of Cbfal mRNA in C3H10T1/ 
2 ceUs.^'' BMP-7 stimulates osteoblast differentiation by increasing col- 
lagen and osteocalcin synthesis, ALP activity and parathyroid hormone 
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(PTH) responsiveness in ROS17/2.8 ceUs and MC3T3-E1 ceUs/""^ 
BMP-7 also increases the formation of bone nodules by rat calvarial 
osteoblasts.''^"^''' Rickard et al. investigated the effects of oestrogen on 
BMP production using two oestrogen-responsive human immortahzed 
osteoblast cell lines, hFOB/ER3 and hFOB/ER9.* Interestingly, oestro- 
gen (17-P estradiol: 10~'''-10~^ mol-L~') increased the expression level 
of BMP-6 mRNA and the production of BMP-6 protein, while BMP-7 
mRNA levels were unchanged. In our study, Bmp7 was upregulated 
more than sevenfold in MC3T3-E1 cells treated with 10~^ mol-L~' 
17-P estradiol for 5 days. The different cell types and culture conditions 
may account for the differences from the previous studies. 

Aside from the role of osteoblasts in bone formation, they also 
express a number of matrix- degrading MMPs, including MMP-2, 
MMP-9, MMP-13 and MMP-14.^''"'^^ Studies oi Mmp null animals 
have defined important roles for osteoblast-derived MMPs. For 
example, Mmp-2 null mice have been shown to exhibit impaired ske- 
letogenesis during development. These unexpected findings can in 
part be explained by the necessity of MMP-2 for osteoblast differenti- 
ation.''^ However, the precise mechanism through which osteoblast- 
derived MMP-2 contributes to osteoblast function is unknown. 
MMP- 14 is an important mediator of MMP-2 activation and could 
promote the secretion of MMP2; both of these are involved in osteo- 
blastic bone formation and/or inhibit osteoclastic bone resorption.^"* 
In this study, the expression of the Mmp-2 and Mmp-14 genes were 
increased in the oestrogen groups compared with the control group. 
This implies that 17-P estradiol likely promotes the proliferation and 
differentiation of MC3T3-E1 cells mediated by MMP-2 and MMP-14. 

CONCLUSIONS 

In conclusion, MC3T3-E1 cells responded to 17-P estradiol depending 
on the concentration and duration of treatment. Based on gene-chip 
analysis, a total of 5 403 genes were determined to be involved in the 
biological activities of osteoblasts either directly or indirectly under 
the influence of 1 7- P estradiol. These findings clearly demonstrate that 
the expression of genes related to osteoblast proliferation, cell diffe- 
rentiation, coUagens and TGF-P-related cytokines increases, while the 
expression of genes related to apoptosis and osteoclast differentiation 
decreases following exposure of MC3T3-E1 cells to minimum essen- 
tial media alpha supplemented with 17-P estradiol. Alternatively, we 
identified several novel genes involved in the response of MC3T3-E1 
cells to 17-P estradiol. Previous studies demonstrated that many of 
these genes were involved in osteogenic pathways in osteoblasts or 
were implicated in the response to 17-P estradiol in other cell types. 

Microarray analysis with functional gene classification plays an 
important role in providing a complete understanding of complemen- 
tary intracellular processes, and the related genes or corresponding 
proteins may be explored in detail to identify more effective methods 
for treating osteoporosis. Furthermore, the obtained biological data 
could be critical for predicting the outcomes of strategies aimed at 
promoting bone regeneration, such as bone tissue engineering. 
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